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Abstract 


Local erosion downstream of hydraulic structures can compromise 
their safety. In the case of works with ski jumps, there are fewer 
empirical studies and empirical formulas to estimate the maximum 
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depth of erosion. Moreover, the uncertainties that still exist in this kind 
of phenomenon, leads to the need to analyze these problems through 
detailed experimental techniques. This paper presents an experimental 
study of local erosion occurs downstream of hydraulic structures low 
altitude with sky jump dissipator in a river gravels and finally the 
experimental results obtained with the empirical formulas of the prior 
art are compared. 


Keywords: Physical model, local erosion, movable bed, Sky jump 
dissipator. 


Resumen 


La erosión local aguas abajo de estructuras hidráulicas puede 
comprometer la seguridad de las mismas. En el caso de obras con 
estructura terminal tipo salto de esquí hay vastos estudios e 
investigaciones disponibles que proponen formulaciones teórico- 
experimentales para estimar la máxima profundidad de erosión. Sin 
embargo, las incertidumbres que aún existen en este tipo de 
fenómenos llevan a la necesidad de abordarlos mediante técnicas 
experimentales detalladas. En este trabajo se presenta un estudio 
experimental de erosión local aguas abajo de estructuras hidráulicas 
de baja altura con salto de esquí, donde la descarga se hace sobre un 
río de gravas; se comparan los resultados experimentales con los de 
fórmulas empíricas que aparecen en el estado del arte, validando el 
grado de aplicabilidad de las mismas. 


Palabras clave: modelo físico, erosión local, fondo móvil, disipador 
tipo salto esquí. 
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The processes of local erosion downstream of hydraulic structures in 
fluvial river beds has been studied by Mason and Arumugam. (1985), 
Lopardo (2005), and Heng, Tingsanchali, and Suetsugi, (2013). They 
arrived at empirical formulations that allow to estimate the maximum 
depth of erosion downstream of hydraulic structures. However, in low 
height structures with ski jumps, there are fewer studies, and 
uncertainties in both the behavior of these hydraulic structures as well 
as in the local erosion process. Generally a solution to particular 
problems is quantified by physical models and the development and 
application of experimental techniques. In this context, physical 
models with mobile beds are used. Where the most difficult parameter 
to reproduce with minimal distortion is the material of the mobile bed. 
For the simplest case, “dragging a particle in a uniform and permanent 
flow” depends on numerous variables such as, the material 
characteristics of the bed, the nature of the fluid, and the mass forces 
acting on the particles (Fuentes, 2002; Lopardo, 2005). 


The objective of this study was to characterize the local erosion 
downstream of low height structures with ski jumps in a physical 
model with a mobile bed. Therefore, it was necessary to perform an 
analysis of the optimum granulometry to represent in the moving bed 
of the physical model. 


Case study 


The case study is Los Molinos derivation Dam, located in the province 
of Jujuy, Argentina (UNC 2012). The dam is located on the river "Rio 
Grande”, upstream of the city of San Salvador de Jujuy and 
approximately 1 km downstream of the confluence with the Reyes 
River (Figure 1). 


The derivation dam is composed of three discharge structures. On the 
left bank of the river there is a fixed weir (Fixed Dam, DF), which is 
separated from the mobile structures by a guide wall. The two mobile 
structures consist of a weir with four circular gates (Mobile Dam, DM) 
and a bottom discharger (or Moderator Channel, CM) regulated by two 
flat gates. Finally, on the right bank of the river is the water catchment 
structure (Figure 2). 
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Moderator Channel 
Guide Wall 


Outlet 
Work 


Figure 2. Aerial view of the structures of Los Molinos dam (Google 
Earth, 2012). 


Since the construction of Los Molinos dam the sediments transport 
have been affected. In the first years, an interruption of the sediment 
flow was generated from the upstream to the downstream sector. 


This process led to the clogging of the reservoir and resulted in a 
generalized descent of the mobile bed downstream of the dam of 
approximately 9 m. This new configuration of the dam environment is 
different from the initial design conditions of the energy dissipation 
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structures. Therefore, the structures began to operate inefficiently and 
local erosions were generated downstream of the structures on the 
right bank of the river. 


To adapt the existing structures to the new configuration, 
modifications were proposed to the dissipation structures which 
consisted in modifying the dissipater type Gandolfo - Cotta of the 
mobile structures by ski jumps. The objective of these structures is to 
move the dissipation zone away from the foundation of the dam so 
that the erosions are produced far from the dam and do not 
compromise the stability of the dam. 


As a complementary structure, to avoid stability problems, a concrete 
wall was designed downstream of the three discharge structures. The 
lower level of the wall was defined as a function of the maximum 
erosion obtained in this study, and the geotechnical characteristics of 
the river material. 


Materials and methods 


Experimental setup 


The tests were carried out in the physical model of Los Molinos Dam 
(Jujuy, Argentina), which was designed to evaluate the hydraulic 
performance of the discharge structures designed to adapt the existing 
structures. The set up consists of a three-dimensional physical model 
(3D) with similarity of Froude, undistorted length scale (E, = 1: 65) 
with fixed margins, and moving bottom in the bed downstream of the 
dam. 


The physical model includes the three discharge structures (Fixed 
Dock, Mobile Dock and Moderator Channel) with their control 
structures (gates). The domain of the model extends upstream of the 
dam, covering the approach sections of the Rio Grande and Reyes 
River with their confluence (approximately 1 000 m) and downstream 
the model extends about 500 m in the prototype. These distances 
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were defined so that the imposed edge conditions do not affect the 
behavior of the flow in the erosion zone. 


To define the geometry of the river in the area of approach to the 
dam, the cross-section method for rivers in regime was applied. This 
method allowed to determine width B, depth H of the cross-section, 
and the slope S as a function of a dominant flow (Farias, 2005). 


To determine the dominant flow, a statistical criterion was applied, 
which is based on the assignment of a recurrence between 1.4 years 
and 2.33 years to this flow. A dominant flow of 600 m?/s was adopted, 
which corresponds to a recurrence of 2 years, and a width of 110 to 
150 m. 


To define the boundary conditions downstream of the dam, a 
numerical model 1D was developed using the HEC-RAS (USACE) 
program. 


Mobile bed material 


In a mobile bed physical model with sediment transport, the results 
depend on the granulometric distribution of the bed material. 


Since it is impossible to achieve the complete similarity of the 
granulometric curve between the model and the prototype, a 
granulometry that best represents the physics of the processes 
involved must be selected to be used in the model (armoring, fluvial 
macro forms, etc.). 


To represent the moving bed in the model, a kind of material with 
density and clastic characteristics similar to those of the prototype was 
used. Figure 3 shows the granulometric curve of the bed at 1:65 scale, 
which will be called the ideal curve. 
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do 0.001 mm | 
dio 0.006 mm 
di 0.009 mm 
da 0.017 mm 
da 0.021 mm 
dao 0.046 mm 
3 0.077 mm 
duo 0.108 mm 
dso 0.215 mm 
de 0.369 mm 
de 0.492 mm 
70 0.586 mm 
d> 0.800 mm 
so 1.000 mm 
des 1.200 mm | | —- scale of model 1:65 
90 1.893 mm — 17 WT | 0 10% 
dos 3.126 mm | | | TTT | ] T T | 0% 
dos 4.689 mm 0.001 0.01 0.1 06 1 10 
dio 9.378 mm . . 
0 11.40 Diameter of the particles (mm) 


Figure 3. Granulometric curve of the prototype downscaled to the 
model ("ideal"). 


To go from the ideal curve to that used in the model, the following 
considerations were taken into account: 


1. The bottom 35 percent of the "ideal" granulometric curve 
corresponds to cohesive sediments (diameters less than 75 microns). 
Therefore, this fraction of the granulometric curve cannot be 
represented in the physical model. 


2. A minimal diameter of the granulometry of 0.6mm is adopted in 
order to avoid the generation of bottom "*parasite” shapes, which 
would modify the bed roughness and would distort the similarity 
between the model and the prototype. 


The selected scale (1:65) allows to represent the upper fraction of the 
"ideal" curve, which includes diameters greater than do. Table 1 
shows that most of the empiric formulas that were developed use 
characteristic diameters included in this fraction of the curve (das, doo 
and dos). 


Table 1. . Parameters of the general formula of the erosion depth 
equation (see citations * in Mason 8 Arumugam, 1985). 
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Selection of the material of the moving bed in the model 


The selection of the material of the moving bed is not trivial and is 
critical for reproducing the phenomenon of local erosion in the 
prototype. 


Three different samples were analyzed, see Figure 4: 


a) Sample 1: Sand with diameters between 1 and 2 mm. 
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b) Sample 2: Sand and Gravel with diameters between 0.6 and 12 
mm. 


c) Sample 3: Sand with diameters between 0.6 and 4mm. 


Scale 1:65 
Sample 1 
=== Sample 2 


— Sample 3 


Diameter of the particles (mm) 


Figure 4. "Ideal" granulometric curve and a curve for each of the 
three samples, UNC (2012). 


In the ideal granulometric curve it is observed that diameters des - dos 
are in the range of 1.2 to 3.4 mm, the three samples analyzed have a 
fraction greater than 25% between these limits. Therefore, it could be 
considered that the three samples may be appropriate to reproduce 
the phenomenon of local erosion. As a consequence of this erosion, it 
is expected that first, that downstream a sediment bar may form with 
the bed material. Second, it is also expected to observe segregation of 
the granulometry without reaching armoring. 


Calibration tests were carried out to determine which of the three 
samples is the one that best reproduces the phenomenology. The 
observations in the three samples were as follows. In sample 1 the 
sedimentation bar was not formed downstream of the erosion zone, 
and erosion continued indefinitely over time. Therefore, the erosive 
process did not stabilize. Consequently, it was considered that this 
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sample did not represent the phenomenon. In sample 2 the sediment 
bar was formed and the erosion stabilized; however, it produced an 
armor in the downstream slope and in the sediment bar with a 
diameter of approximately 8 to 12 mm, equivalent to 520 and 780 mm 
in the prototype. The material of this size in the prototype is a small 
proportion (less than 5%). This is the reason why it was considered 
that this type of armoring could not be formed. Thus, sample 2 was 
rejected since it was not representative of the phenomenon. In sample 
3 the sediment bar was formed and the erosion stabilized in addition, 
there was a segregation of the sediments in the downstream slope, 
and in the sediment bar with an average diameter of 4 mm, equivalent 
to 260 mm in the prototype (Figure 5). It was considered that this 
sample was the one that best reproduced the physical phenomenon, 
and was the one chosen to carry out the erosion tests in the model. 


Figure 5. Sediment bar downstream of the erosion zone, Sample 3- 
Qprototype = 900 m/s. 


It is concluded that, regarding the case analyzed, the stabilization of 
the erosion depends on the formation of the sediment bar, and it 
depends on the transport of sediments. Therefore, the dimensions and 
shape of the local erosion is a function of the hydraulic and sediment 
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characteristics as well as their granulometric distribution. Interestingly, 
in this work the final erosion is a function of the granulometric curve 
unlike Mason and Arumugam (1985) and Heng et al. (2013) where the 
sediment bar formation is not considered. 


Modeled hydraulic conditions 


Test 1. Flow = 900 (m?*/s) spilling through the moderator channel and 
the mobile dam. 


Test 2. Flow=3200 (m*/s) spilling through the fixed dam. 


Test 3. Flow=4200 (m*/s) spilling through the moderator channel, the 
mobile dam, and the fixed dam. 


Test 4. Flow=1600 (m/s) spilling through the moderator channel, the 
mobile dam, and the fixed dam. 


Test 5. Flow=600 (m*/s) spilling through the moderator channel and 
the mobile dam. 


Experimental methodology 


The erosion zone was surveyed with an optical level in three profiles 
parallel to the structure but perpendicular to the flow. Profile 1 was 
measured next to the dam, profile two was taken at the maximum 
depth of the erosion pit, and profile three was obtained on the 
sediment bar. 


Also, the erosion surface was measured (Figure 6a) using a digital 
camera (Microsoft, 2010) which allows the acquisition of planimetric 
information with high spatial resolution. With this information a digital 
elevation model was generated (DEM, Figure 6b) similar to Pagot et al. 
(2014). 
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Canal moderador y dique móvil ¿¿ 


Figure 6. (a) Optical image of the erosion downstream of the mobile 
dam using Kinect, y (b) Digital Elevation Model. 


This survey technique has the advantage of acquiring high spatial 
resolution data remotely without being in contact with the surface. 


To ensure that the erosion pit reaches a state of dynamic quasi- 
equilibrium, erosion monitoring was carried out while the tests at 12 
control points were carried out. During the first 90 minutes of the test, 
measurements were taken every 15 minutes and then every 30 
minutes until the erosion was stabilized. It was considered that when 
three consecutive measurements did not differ significantly, this 
equilibrium state was reached. 


In addition, to locate the zones of maximum velocity and the 
transversal distribution of the flow, the entrance velocities to the 
structures were measured by PTV technique (Velocimetry by particle 
tracking). The velocimetry technique by particle tracking allows to 
obtain the instantaneous and average velocity fields, trajectory lines, 
and average trajectory lines. To analyze the images, the PTVLab 
software (Patalano, 2017) was used. 


Empirical Estimates of the Maximum Depth of Erosion 


(hs) 


There are numerous investigations that present semi-empirical 
formulations that allow to determine the maximum depth of erosion 
based on energetic considerations of the flow and the characteristics of 
the material of the bed such as typical diameters. 
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Below are the equations used in this work to determine the maximum 
depth of erosion. These equations were classified into three groups 
according to the parameters that intervene in each of them: 


Group I -— Equations that express the maximum depth of erosion (h,) 
in terms of the energy difference between upstream and downstream 
(H), the specific flow of the jet (q) and the characteristic diameter of 
the bed material (d). 


The general formula of this group of equations is as follows: 


qua? 
dz 


where: K, x, y, z are coefficients and exponent that were calibrated by 
different authors, in Table 1 the values of the coefficients proposed by 
each author are presented. 


Formulas 1 to 11 have zero exponent of the characteristic diameter. 
Therefore, these expressions do not take into account the diameter of 
the sediments to determine the maximum depth of erosion pit. 


Within this group the formula proposed by the 1.D.I.H is included: 
Formula 23 - 1.D.I.H. (1990): h, = 2.6662qH705 + 0.3291H 


Group II - These equations contemplate in addition to the variables q, 
H y d, the height of the water at the point of impact (h,.). 


h, 0.333 
Formula 24- Jaeger (1939): h, =0.6q0504J025 (52) 


Formula 25 - Martins — A (1973): h,=0.14N — 0.7324 17h, 


Q3H15 
a? 


where: N=7 


The discharge Q can be substituted for q because it is assumed that 
this expression may be applicable for large flow sheets. 


0.60 y0.05 0.15 
Y 


Formula 26 - Mason (1984): h, = 1.652 


a0o-10 
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Group III - These equations also consider the angle at which the jet 
hits the water (B). 


_ % ] _ 1.804qsinf 1 _ 1126 
Formula 27 - Mikhalev (1960): Ms on (BP 2) 
Formula 28- Rubinstein (1963): 

0.75 1.20 
hs = hop + 019 (z) (csm) 

_ s _ _ [0.97 _ 135 qsinf 

Formula 29 - Mirskhulava - (1967): is = (E A 


0.25h, 


The 29 formulas considered are semi-empirical and were calibrated 
from experimental data obtained from physical models or erosion 
values observed in the prototype. 


To estimate the specific flows (q) that each of the structures spill, two 
aspects were considered. First, the H-Q curves corresponding to each 
structure, and second the flow was distributed uniformly along the 
width of the structure. The height of energy (H) was obtained by 
adding the geometric difference between the level upstream and the 
restitution to the speed energy in the section in which the 
measurement of the levels was carried out. 


2 
(2x9) 
where, 


H¿: geometrical difference between the upstream level and the level of 
restitution, 


V, : average speed in the upstream section. 


To assign values to the characteristic diameters that intervene in the 
formulas the granulometric curve of the bed was used (dos=0.130m; 
do0=0.123m; da5=0.116m y d:=0.068m). 
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Results 


Experimental studies 


The maximum depth of erosion (h,) was plotted against the two 
variables that appear in all the empirical formulations: the specific flow 


q" and the height of energy "H". 


In principle, it is expected that as the specific flow rate (q) and the 
power height (H) increase, the maximum depth of erosion (h;,) will also 
increase (Heng et a/l., 2013; Ahmed-Amin, 2015). However, in Figures 
9 and 10 where the relationships q -h, and H-—h, are plotted, a 
different behavior was observed for the Moderator Channel and the 
Mobile Dam. Figure 9 shows that in the Mobile Dam the maximum 
erosions (h,) did not present significant changes between the flows of 
13.4 m?/s and 16.2 m/s. In contrast to the expected increase in (h,) , 
in the Moderator Channel an erosion of 15 m was observed for a flow 
of 18 m/s, and an erosion of 14m for a flow of 19.7 m/s. Similarly, in 
Figure 10 there were no variations in maximum erosion when H 
increased from 12.4 m to 13.6 m. 


To explain the lack of monotonous behavior, two hypotheses were 
formulated: 


Hypothesis1: The lack of monotonous behavior is produced by a 
change in the flow direction upstream of the structures. In the tests 
carried out in the physical model, a change in the flow direction was 
observed when the water drained by the Fixed Dam. 


Hipótesis2: The behavior observed in Figure 9 and Figure 10 is due to 
a change in the H-Q relations of the Moderator Channel and the Mobile 
Dam. When the water level exceeds the height of the gates they begin 
to function as an orifice. This is the reason why there is a change in 
the H-Q relations. 
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Approach flow 


Applying the PTV technique (Patalano, 2017) it was observed that, 
when the water drains only by the Mobile Dam and the Moderator 
Channel, the flow current lines enters obliquely (Figure 7a). When the 
water drains through the three structures, the current lines upstream 
of the dam enter perpendicular to the dam (Figure 7b). 


Figure 7. Mean trajectory lines upstream of the structures: (a) 
discharging only by the Mobile Dam and Moderator Channel, and (b) 
discharging through the tree structures. 


Discharge curves 


To quantify the discharges of the structures, their height-flow curves 
(H-Q) were obtained. The fitting was made with experimental data 
obtained from the physical model for each of the structures. Figure 8 
(a) shows the height-discharge curve of the Mobile Dam and in Figure 
8 (b) that of the Moderator Channel. It can be seen that the Mobile 
Dam begins to work as an orifice when the height above the landfill ¡is 
approximately 5.4 m and the Moderator Channel for a height of 6 m. 
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Figure 8. Height-discharge flow curves (H-Q) of the Mobile Dam (a) 
and Channel Moderator (b). 


Figure 9 shows, in a dashed line, the specific flow corresponding to 
Hypothesis 1 and in a continuous line the flow corresponding to 
Hypothesis 2. It is observed that the change in the flow direction 
upstream of the structures (Hypothesis 1) has no significant influence 
on the unexpected behavior of the relationships (q-hs) and that the 
specific flow corresponding to Hypothesis 2 (change in the height- 
discharge relationships), for the two structures, is in the range of flow 
rates in which a different behavior than expected was observed. Then, 
it is concluded that the unexpected behavior between the discharge vs. 
erosion depth relationships can be explained by the change in the H-Q 
relationships. 
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Figure 9. Experimental relationship q-hs of CM, DM and DF. 


In Figure 10 it is observed that the change in slope in the relationships 
h¿-H occurs for H values between 12.4 m and 13.6 m for both 
structures. The value of H for Hypothesis 1 is 12.3 m for both 
discharge structures and 11.8 m, and 13 m for Hypothesis 2. 
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Figure 10. Experimental relationship H-hs of CM, DM and DF. 


Empirical estimates 


The erosion values (h,) obtained by applying the 29 empirical formulas 
showed a great dispersion in the results. The standard deviation value 
is approximately 50% of the average value of the 29 formulas. 


Next, the experimental results were compared with the estimates. 
Figure 11 (a), (b) and (c) present the scatter diagrams of the 
experimental erosion depth versus estimated erosion for the Moderator 
Chamnel, the Mobile Dam, and the Fixed Dam, respectively. It was 
observed that no formula worked satisfactorily for the three structures. 
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Figure 11. Comparison diagram between experimental erosions and 
empirical estimates, (a) Moderator Channel, (b) Mobile Dam and (c) 
Fixed Dam. 


Figure 12 shows a BOX-PLOT graph in which the results for the three 
structures are presented. It can be seen that the experimental result 
for the Moderator Channel is close to the second quartile (Q2) and that 
the experimental values corresponding to the Mobile Dam and the 
Fixed Dam are located close to the third quartile (Q3). It is concluded 
that the maximum depth of erosion not only varies with the 
sedimentological conditions, which in this case remain constant for the 
three structures, but also with the hydraulics of each work in 
particular. 
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Figure 12. BOX-PLOT graph of Chanel Moderator, Mobile Dam and 
Fixed Dam for a flow rate of 4 200 m/s. 


Table 2 shows the percentage relative differences between the 
maximum erosions measured in the physical model and those 
estimated with the different empirical formulations for each of the 
structures and each of the scenarios tested. The formulas that were 
best adapted to the conditions tested for the three discharge 
structures were those proposed by Zeller (1981) and Machado B 
(1982). Those that best adapted to the experimental measurements 
made in the Moderator Channel were Rubinstein (1963), Lopardo 
(1987), Mirskhulava (1967), Zeller (1981), Machado B (1982) and 
INCYTH (1985); for the Mobile Dam were Veronese (1937), Sofrelec 
(1980), Machado (1982), Mason and Arumugam (1985), IDIH (1990), 
Jaeger (1939) and Zeller (1981) and for the Fixed Dam were Machado 
B (1982), Veronese (1937), Zeller (1981), Mason and Arumugam 
(1985) and Sofrelec (1980). 


Table 2. Percentage difference between the empirical estimates and 
experimental erosion for the CM, DM and DF. 


Moderator channel Mobile dam Fixed dam 

Discharge total 160 | 420 | 60 160 |420| 160 | 420 | 320 
(m/s) E cis 0 0 0 200 0 0 0 0 0 
Esien 10. | 15. |14.3|15.4[9.9| 13. [13.8 17.7|10.4|14.8 | 18.9 
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Experimental 5 2 9 
hs (m) 
Percentage difference between the empirical estimates and 
Formula ; Ñ 
experimental erosion (%) 
Veronese 
1 (1937) 44 5 19 22 -3 -2 12 1 12 -4 -15 
2| PESA as .6 | 26 | 23 [230139 300) 36 08: |<37 [44 
(1966) 
3] Pame-B [os laz | <38 |:<36 |=49 250.|-42 | 246: | <40: |:<48 1 253 
(1966) 
a Pame e ls caiono. sz (66.66 | 61 ¡64 | 60: |-<65 | 63 
(1966) 
Chian Min Wu 
5 (1973) -16 |-39| -30 | -29 |-42| -42 -35 | -41 | -33 | -44 | -50 
Taraimovich 
6 (1978) -26 |-46 | -37 | -35 |[-55| -51. -43 | -47 | -47 | -52 | -57 
Machado - A 
7 (1982) 114 | 56 77 81 50 | 47 67 50 71 45 27 
8 | Sofrelec (1980) | 52 | 10 | 25 27 | -2 1 15 5 12 -4 | -15 
9 | INCYTH (1985) | 2 |-26| -16 | -14 |-29| -30| -21 | -29 | -19 | -31 | -40 
1 Martins - B 
0 (1975) -1 |-28| -19 | -17 |-36|-34 | -25 | -32 | -27 | -37 | -44 
1 | Lopardo et al. _ , FOR - - - - - - 
a (1987) 6 22 9 6 26|-26 | -14 | -21 12 23 32 
1 Schoklitsch 
> (1932) -21 | -42 | -34 | -33 |-48| -47 -39 | -45 | -40 | -48 | -54 
1 Veronese 
3 (1937) -53 |-66 | -61 | -60 |-68| -68  -63 | -67 | -63 | -69 | -72 
1 Egenberger ' - _ ' , Ñ Ñ E S - q 
4 (1943) 7 31 | -19 15 |-40|-36 | -25 30 27 34 | -40 
: Hartung (1959) | -11 | -34 | -23 | -20 1-44 | -40|-30 | -34 | -33 | -40 | -45 
A Franke (1960) | -31 | -48 | -39 | -35 |-58| -53 | -44 | -47 | -48 | -52 | -56 
, Kotoulas (1967) | 191 | 117 | 153 | 167 | 73 | 93 | 129 | 118 | 111 | 93 78 
a Zeller (1981) 20 | -9 9 18 |-28|-18 | O -3 -8 | -13 | -18 
1 | Chee y Padiyar 
9 (1969) 147 | 81 [108 |116|50| 63 89 | 76 | 78 57 | 42 
2 |Bisaz y Tschopp 
0 (1972) 90 |39| 58 | 61 |30| 30 | 48 | 34 50 | 28 13 
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2 Chee y Kung 

1 (1971) 83 | 34 53 58 [18 | 23 | 41 30 38 19 7 

2 Machado - B 

> (1982) 30 | -5 9 11 | -9 | -10 3 -7 5 -10 | -21 

E 1.D.1.H (1990) 61 17 31 35 -4 2 17 9 10 -5 -15 

A Jaeger (1939) | 43 7 22 | 44 0 1 15 20 17 15 -4 

2 Martins - A 

: (1973) -28 | -45 | -35 | -13 |[-54|-51 | -41 | -29 | -44 | -35 | -47 

> Mason y 

6 Arumugan 45 | 6 |19|30|-7|-3|10| 7 7 -2 |-15 
(1985) 

- Mikalev (1960) | 47 9 27 37 |(-57|-43| -31 | -31 | -61 | -63 | -65 

A A Flats 1 [ss sz | .20 | 02 [ss 32 (edo 

2 Mirskhulava 

9 (1967) 20 |-10 5 13 |-87|-83 | -79 34 49 63 62 


Finally, the experimental results observed in this study are compared 
with those presented in similar cases using the graph by Lopardo 
(2005). In this graph, the main variables are related in a 
nondimensional way. Figure 13 shows the curves corresponding to the 
expression proposed by INCYTH, together with experimental points of 
erosions in prototype works (Ptos exp.), the envelope curve given by 
the design expression proposed by Lopardo, and the results obtained 
in this work for each of the structures. It is observed that the points 
measured in this work are close to the envelope proposed by Lopardo 
and in some cases h,/H is higher. 
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Figure 13. Non-dimensional graph of the variables considered. 


Conclusions 


Tecnología y ciencias del agua, 10(3), 70-95. DO!: 10.24850/j-tyca-2019-03-04 


0.4 


93 


¿ 2019, Instituto Mexicano de Tecnología del Agua 
Open Access bajo la licencia CC BY-NC-SA 4.0 
Tecnología y (https: //creativecommons.org/licenses/by-nc- 


CienciaszAgua sa/4.0/) 


In the calibration stage "Selection of the granular material of the 
mobile bed" it was observed that the erosions depend on the formation 
of the sediment bar downstream of the erosion pit and the formation 
of the bar depends on the characteristic diameters of both the particles 
and the granulometric curve. 


When the jet of water hits the bed, a great dissipation of macro- 
turbulent energy is generated, which initiates the local erosion 
process. The bottom material is removed, placed in suspension, and 
transported downstream where the flow reduces its energy, and 
therefore its transport capacity. 


Finally, the sedimentation of the larger particles takes place, and 
initiates the formation of the sediment bar. The latter contributes to 
stabilize the erosions, reaching lower depths than those that would 
occur without the formation of the sediment bar, due to two effects. 
First, it increases the level of restitution in the energy dissipation zone. 
This reduces the value of H and increases the height of the water in 
the impact zone h,, Second, it raises the level of the bed downstream 
of the erosion pit. Therefore, for the sediments to come out of the 
erosion pit they must be raised to a larger height. 


Regarding the "analysis of experimental results", it is concluded that, 
for structures of low height regulated by gates, the relations between 
q and H with the maximum depth of erosion (hs) vary with the 
discharge law (H-Q) of the structures. A large dispersion was observed 
in the results when erosions (hs) were calculated with "empirical 
formulas". The application of these formulas must be limited to 
preliminary estimates to obtain an order of magnitude of the potential 
erosions. 


When only empirical expressions are used, formulas that have been 
calibrated for hydraulic and sedimentological conditions similar to 
those used in this work should be applied. In this study, the formulas 
that showed the best behavior for the conditions tested are those by 
Zeller (1981) and Machado-(1982). 
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